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The Lax representation for the nonstationary Schrodinger equation with rather arbitrary 

potential is proposed. Some examples of the construction of exact solutions are given by means 

' of Darboux Transformation method. 
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1. Introduction. 

The increase of interest to the non-stationary Shrodinger equation (NtSH) noticed dur- 
ing last decades is apparently explained not only by needs to solve quantum-mechanical 
problems, but also by its role in the method of inverse problem (for example, in the theory 
of the Kadomtsev-Petviashvili equation) and related problems [1-4]. 

It is necessary to note, that known methods of construction of the exact solutions of 
NtSH equation are basically concentrated around two methods: the method of intertwin- 
ing relations (its special case is a method of Darboux Transformation (DT)) [1,2] and 
the method of dressing chains (see, for example, [3]). These methods have both some 
advantages and definite defects (without details, we shall specify only, that they do not 
give a regular procedure of construction of the solutions neither for any rather "good" 



O ■ potentials, nor for the case of higher spatial dimensionality). 



It is necessary to mention separately the paper [5], in which the scattering theory was 
developed directly for the NtSH equation, and its connection with the solutions of the 
Kadomtsev-Petviashvili equation was discussed. 
Ch ■ In a certain sense the alternative approach was proposed in [6]. It is based on the 

> . application of the Darboux-like anzatz, it does not depend on the spatial derivatives (and, 

^ \ hence, on the dimensionality of the problem), and it also reduces the solution of NtSH 

^ \ equation to the solution of some functional-differential equation, which does not contain 

any potential. 

The idea of the present paper is the following. The NtSH equation (though it is 
linear, but it includes a coefficient function - potential) is considered as a (rather special) 
particular case of the nonlinear Shrodinger equation: 

= + x(^, ^)^, (1.1) 

assuming that x(\l/, ^f) = W{x^ t), where W{x^^) is a given function depending only 
on spatial and time arguments. Then this equation represents the completely integrable 
system. It means, that it possesses also all appropriate attributes of such systems: the 
scattering problem, an infinite and involutive series "of the conservation laws" (since 
the initial system is not conservative, no sense to speak about "the relations of balance"), 
identities of traces, hierarchy of the higher equations etc. Certainly, besides that it also 



1 



allows to construct the exact solutions of the initial equation, for which only limited class 
of solutions is known up to now (see, for example, monography [7]; let's also note here 
the important review [8], in which the quantum-mechanical systems were discussed from 
the point of view of the representation theory of groups and algebras). 

Thus, the problem is reduced to the construction of a Lax pair and to the analysis of 
the appropriate associate linear system. 

The paper is organized as follows. In Section 2 we obtain the Lax representation by 
means of transition to the Heisenberg picture. Section 3 is devoted to verification of the 
covariance property and to calculation of the exact solutions by the DT method. 

2. Lax Representation and Heisenberg Picture. 

Let us consider the NtSH equation (here and below we use the system of units, in 
which the Plank's constant h — 1, and the mass of particle m = 1/2) in the coordinate 
representation: 

i^t = -'^.. + W^, (2.1) 

where \E' = \E'(a;, t) is the wave function ^ Potential W = W{x,t) is supposed to be given 
(it plays a role of functional parameter) and real, therefore: 

im%^i^,^-^,^%. (2.2) 

Let's rewrite the equation (2.1) as a condition of compatibility of linear system of (2 x 2) 
matrix equations: 

$^ = t/$A, $i = V$A, (2.3) 

where $ = $(,t, t, A) e Mat{2, C), A = diag{X, A), A e C is the parameter, ImA ^ 0, and 
let's assume that 

Here Uq = Uo{x,t), Vq = Vo{x,t), fo = /o(a;,t), go = go{x,t), po ^ po{x,t), qo = qo{x,t), 
and the equation (2.1) is equivalent to the relation: 

uot = ^'ox- (2.5) 
It can be considered as the constraint, which can be solved as: 

a /■* 

Uo = i'^ix,t) / W{x,t')'^{x,t')dt' + Ao{x), 

2 Jo 



(2.6) 



2 - a r 

vo = -^x(x, t) + — — / W{x', t)^{x', t)dx' + Bo{t), 



^As was noticed in [5], the replacement ^ 4'cxp(At/i) reduces the equation i'^t = —'^xx + 
{W{x,t) — A)\E', where A G C, to the form in (2.1), and thus the analysis given below is applicable also 
to this equation. 
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where a G M is the parameter, Aq{x), BQ{t) are arbitrary functions, and M > is 
arbitrary (finite or infinite) number, introduced in order to provide the convergence of 
integrals below. 

The requirement of compatibility of system (2.3) leads to relations: Ut — Vx — and 
[U, V] = 0; by the straightforward calculations one can check that they will be fulfilled, 
if the functions /o, go will satisfy the linear equations in partial derivatives of the first 
order: 



(ln/o)t (ln/o)j; = (ln^/;o)x, (Ing'o)t (Ing'o)x = (Inwo)a;, (2.7) 

Wo Wo Wo Wo 

where wo = (uq - uq)/{vo - wq = wq. 

For the solution of these equations we shall introduce the real function as the 
curvilinear integral: 

6 = t) =i {uo- Uo)dx + {vo - Vo)dt, = ^i, (2.8) 

{xo,to) 

where xq is the initial point, to is the initial moment of time (later on we put to = 0). 

By virtue of a condition (2.5) the integral (2.8) does not depend on the integration 
path. It allows to integrate the equation (2.7). Taking into account that Wo = Cix/Cit, 
we have: 

/o = 6.^o(ei), ^0 = 6x^0(6), (2.9) 

where Fo{.), Go{-) are arbitrary functions. Therefore po — CuFoiCi), Qo — CitGoiCi)^ 
finally we obtain: 

(2.10) 

Thus, using the Heisenberg picture, we formally have found the Lax representation 
for equation (2.1). Let us notice that the requirement of reality of potential together 
with the condition (2.2) does not impose any additional restrictions on matrix elements 
of matrices U and V. 

Let us consider in more detail the Heisenberg picture, which appeared here. For that 
the equation (2.1) can be rewritten as: 

i^f ) = H{t)¥^\ (2.11) 

where H{t) is the Hamiltonian, assumed hermitian and time-dependent. Let's assume 
that 

^(5) = U{t)¥^\ (2.12) 

where xj/W are the vectors of state in the pictures of Shrodinger and Heisenberg, 

correspondingly, U is the unitary operator. From (2.11), (2.12) it follows, that this 
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operator satisfies the equation: tjt = —iH{t)U, whence by the assumption, that U{0) — I, 
where / is the unit operator, we obtain: 



U{t)^fe-'^oH{t')dt'^ (2.13) 

where T is the time-ordering operator (Dayson operator). 

According to the basic principles of the Quantum Mechanics, if L°, L{t) are the 
observables in pictures of Shrodinger and Heisenberg correspondingly, the average values 
of these observables in both pictures coincide: 

< ^(5)|i,0|^(5) *W|L(t)|^(^) > . 

Then taking into account the unitarity of U, from the last equation we shall obtain 

L(t) = U+{t)L^U = IJ-^L^U. (2.14) 

where the symbol " + " means the hermitian conjugation. Differentiating (2.14) in time, 
we shall have: 

Lt = i[U-^HU, L]. (2.15) 

It is possible to obtain the representation equivalent to this, using the another, in com- 
parison with (2.13), relation for the operator U. Actually, it is easy to see, that 

U{t) f H{t')U{t')dt' = 

Jo 

°° 1 /1\ " /■* /"* /■* 
= V -r - dt, dt2... dtnT{H{h) . . . H{tn)). (2.16) 

n=0 \^/ "^0 Jo Jo 

Let us assume that R — [H, U]. Then due to (2.16) we have: 

°° 1 /"* /■* /■* 

- dh dt2... dtn[H{t),f{H{h)...H{Q)], 

n=0 ^ ^ ^ J^ J^ J^ 

and, hence, 

Lt^i[H + U-^R, L]. (2.17) 

Thus, equations (2.15) and (2.17) are the generalizations of Heisenberg's equations of 
motion for the case when operators H{t) and U do not commute; and the equation (2.17) 
is equivalent to (2.1). It is easy to show, supposing that the quantity Uq in (2.6) is the 
observable in a picture of Heisenberg, setting formally uq = Ki'^^^\ vq = K2'^^^\ where 
Ki, K2 are the integral and integro-differential operators correspondingly, and obtaining 
the closed equation for Uq: Uot — [K2{Ki^Uo)]^. 

3. Darboux Transformation and Exact Solutions. 

Let us proceed to the construction of exact solutions of equation (2.1). With that end 
in view, we shall use a method of matrix Darboux Transformation [9], setting in (2.6) 
a = 0. 
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Let $1 is the solution of system (2.3), corresponding to a choice of solution ^ = of 
equation (2.1) with some fixed value A = Ai. Let us assume that Li — $iAj"^$j"\ Li e 
Mat{2, C). Then we shall check the covariance of the system (2.3) under the transforma- 
tion: 

$ - $ - Li$A. (3.1) 
This procedure results to the following "dressing" formulas: 

u = uuL^\ [/ = t/ + $i[Ar\$r'$ix]$r\ 

(3.2) 

V = L^vL^\ v^v+ $i[Ar\ $r^$u]$r\ 

each of which, as it can be checked by direct calculations, is equivalent to others. 

The formulas (3.2) can be slightly simplified, using possible symmetries (involutions) 
of the problem. Actually, combining the relations A = aiKcxi and A = o'2Ao'2, where 
(Ti, cr2 are Pauli matrices, and correspondingly: U = aiUai (if po = Qo) and U = (72^7(72 
(if Po — —Qo), we obtain $ = (Ti$(Ti, and $ = (72$(72, correspondingly. 

From the second of relations (3.2) after some calculations we have: 

^^^ + iY,, y = Y(x, t, AO = (r^l^) ' ^^-^^ 

where y = y{x,t,Xi) = (p^^ {x,t, Xi) / (p^^i {x,t, Xi) , (p\^\x,t, Xi) are matrix elements of the 
matrix $i. 

The expression for "dressed" potential of the NtSH equation can be found either by 
replacing in (2.1) \i/ — > ^ and W W, or by using the following reasons. From fourth 
of relations (3.2), analogously (3.3), we have: 

Vo = Vo- Yf (3.4) 
From the explicit expression for Vq we find: 

*'(M)^mM)^-^->^. (3.5) 

Let us notice also that the reality of potential W follows from (3.3). 

Now for construction of the solutions of equation (2.1) it is necessary to integrate the 
system (2.3). For this purpose we shall rewrite it as the overdetermined scalar system of 
the Riccatti equations for the function y: 

Vx = -Ai[/o2/^ + {uo - uo)y - /o], yt = -Ai[poy^ + {vo - vo)y - po]- (3.6) 

Therefore in particular it follows, that y satisfies the linear equation in partial derivations: 

yx = woyt, which has the solution of a kind y = R{^i), where /?(.) is an arbitrary complex- 
valued function of ^i, determined by relation (2.8). To find this function, we shall consider 
first of the equations (3.6) at t = ^. For its solutions we shall suppose, that y = at 

^As it follows from (2.8), the corresponding expression is ^i(a;,0) = i J^j^{uo{x' ,0) — uo{x' ,0))dx. 
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X = —M and y is a bounded function for all x > M . Let us take advantage of an 
arbitrariness in a choice of function /q. Then in the simplest case of integrability it is 
necessary to choose /o = —Jo = Cq = Cq = const. Coming back to the variable for the 
function y, we shall find ( Ai = — Ai, y{x, t) — y{x, t)): 

Y{x,t) = = J^"",^' , ei = ei(6) = r [uo{x',0)-uo{x\Q)]dx'. (3.7) 

v2 + tan©i V 2 7-M 

Thus, from (3.3), taking into account (3.7), we shall have: 

^ ^^ + 2i["°^^' 0)-^o(^, 0)]No(^, t)-Mx, ^)](v/2 + tanhe0. (3.8) 
cosh2ei(\/2tanhei + l)2 

The Darboux Transformation (3.1) can be multiply iterated, leading multiply "dressed" 
solutions (procedure of their construction for solitons of deformed Heisenberg magnet 
based on the dressing formulas similar to (3.2) is given, for example, in [10]). 

Let us discuss two simplest examples, illustrating an application of the proposed ap- 
proach. 

a) . Let W{x,t) = 0, i.e. we shall consider the free NtSH equation, and let us choose 
^ = as a plane wave: = ae^'^^^^^ *, where a is the amplitude of the wave, k is the 
wave number. Then, taking for simplicity M = 0, from (2.8) we have: 

4a 

6 = — -sm{kx - kh), (3.9) 

rv 

and the expression for "dressed" wave function reads: 

cosh2ei(V2tanhei + l)2 

where 

6i = ©i(6) = 

This wave function will correspond to the potential W{x, t) — W{^i), determined from 
(3.5). Let us notice here an occurrence of the harmonic with the "opposite" phase in the 
spectrum of excitations of the system. 

b) . Let the initial solution be the well-known wave function for the harmonic oscillator 
( W{x) = {1/A)ujIx^) (see, for example, [11]): 

^U^,t) = c.e-^-("+V2)*e-^i7„(y^,;), (3.11) 

where uq is the frequency, -ff„(.) is Hermit polynomial, c„ are the real normalization 
constants. Taking ^'i„(a;, t) = X{x)T{t) and X = X, according to (2.8), we shall find 
(M = oo): 

/X pt pt px 

Xdx-iX,l,=-oo / {T-f)dt + i / {T-f)dt / WXdx, (3.12) 
-oo Jo Jo J-oo 
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and 



©1 = = -^c„Ai / e 4 HM —x)dx, Ai = -Ai. 

-M V ^ 



Then, from (3.3) we find 



^2c„i/'„(y^x)e-i'^o^' cos[o;o(n + i)t](V2 + tanhOi)^ 



(3.13) 



cosh^ei(V2tanhei + l)2 

and, hence, we have here the same effect, as in the previous case. 

In conclusion of this section we shall consider the connection of Darboux Transforma- 
tion and the theory of the coherent states. For this purpose we shall represent the wave 
function of this state, which is defined as the minimizing for the uncertainty relations ^, 
as a decomposition in the wave functions of the harmonic oscillator: 

*cofeer = Xl^n*ln, (3.14) 

n 

where '^in are defined by the expression (3.11), and are coefficients of the decomposi- 
tion. After the single application of Darboux Transformation we have: A^ — > A„, — > 
and, besides that, 

= (3.15) 

n n 

where 

= y" '^coher{x,t)^ln{x,t) dx, ^ J coher{x,t)^ln{x,t) dx. 

Then taking into account (3.3), we obtain the connection of coefficients A^ and A^. 

An = An-i J ^'coher^'a; dx, (3.16) 

where the quantity Y is determined in (3.3). 



4. Conclusion. 

Let us notice, that it is possible also to find a slightly different Lax representation 
for the equation (2.1), allowing to study the scattering theory for the operator of the 

^It is possible to define the coherent states in another way as eigenstates of non-hermitian annihilation 
operator; then the quantities A„ are connected with matrix elements of this operator. 
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associated linear problem. On this basis it is possible to obtain, in particular, "the re- 
lations of balance" (instead of conservation laws, because an initial system, as it was 
already marked, is not conservative), and also to construct the hierarchy of the "highest" 
equations. It is especially interesting, that already the simplest equations appear to be 
nonlinear. 

As a whole, such approach can probably serve as a basis for quasiclassical quantization 
of NtSH equation as a completely integrable system. We are going to present these results 
later on. 

The author is grateful to I.V.Komarov and M.V. loffe for the useful critical remarks. 
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